Overpopulation of in pp collisions: 
a way to distinguish statistical hadronization from string dynamics 
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The fi/H ratio originating from string decays is pre- 
dicted to be larger than unity in proton proton interactions 
at SPS energies (_Ei a b=160 GeV). The anti-omega dominance 
increases with decreasing beam energy. This surprising be- 
havior is caused by the combinatorics of quark-antiquark pro- 
duction in small and low-mass strings. Since this behavior is 
not found in a statistical description of hadron production in 
proton proton collisions, it may serve as a key observable to 
probe the hadronization mechanism in such collisions. 

Hadron yields and their ratios stemming from the final 
state of ultra-relativistic heavy-ion collisions have been 
extensively used to explore the degree of chemical equi- 
librium jl]-[)) and to search for evidence for exotic states 
and phase transitions in such collisions jjj . Under the as- 
sumption of thermal and chemical equilibrium, fits with 
a statistical (thermal) model have been used to extract 
bulk properties of hot and dense matter, e.g. the temper- 
ature and chemical potential at which chemical freeze-out 
occurs [||-[|. 

The application of a statistical model to elementary 
hadron hadron reactions was first proposed by Hagedorn 
in order to describe the exponential shape of the 
TOt-spectra of produced particles in p+p collisions. Re- 
cent analysis pd| ] on hadron yields in electron positron 
and proton proton interactions at several center-of-mass 
energies have shown that particle abundances can also 
be described by a statistical ensemble with maximized 
entropy. In fact, they are consistent with a model as- 
suming the existence of equilibrated fireballs at a tem- 
perature T «160-170 MeV. These findings have given 
renewed rise to the interpretation that hadronization in 
elementary hadron hadron (hh) collisions is a statistical 
process, which is difficult to reconcile with the popular 
picture that hadron production in hh collisions is due 
to the decay of color flux tubes a model which has 
explained many dynamical features of these collisions. 

In this letter we argue that the O/O = + /0~ ratio 
in elementary proton proton collisions is an unambigu- 
ous and sensitive probe to distinguish particle produc- 
tion via the break up of a color flux tube from statistical 
hadronization. 



Color flux tubes, called strings, connect two SU(3) 
color charges [ 3 ] and [ 3 ] with a linear confining poten- 
tial. If the excitation energy of the string is high enough 
it is allowed to decay via the Schwinger mechanism [|13[ , 
i.e. the rate of newly produced quarks is given by: 



dp± 
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where k is the string tension and m_i_ = \Zp'± + m2 is the 
transverse mass of the produced quark with mass m. 

However, specific string models may differ in their phi- 
losophy and the types of strings that are created: 

• In UrQMD juj the projectile and target pro- 
tons become excited objects due to the momentum 
transfer in the interaction. The resulting strings, 
with at most two strings being formed, are of the 
diquark quark type. 

• In NeXuS M|, the pp interaction is described in 
terms of pomeron exchanges or ladder diagrams. 
Both hard and soft interactions take place in par- 
allel. Energy is shared equally between all cut 
pomerons and the remnants. The endpoints of the 
cut pomerons (i.e. the endpoints of the strings) 
may be valence quarks, sea quarks or antiquarks. 

• In PYTHIA a scheme similar to that in 
UrQMD is employed. However, hard interactions 
may create additional strings from scattered glu- 
ons and sea quarks. Most strings are also of diquark 
quark form. 

Fig. [I] (left) depicts the anti baryon to baryon ra- 
tio at midrapidity in proton proton interactions at 
160 GeV. The results of the calculations by NeXuS, 
UrQMD and PYTHIA, which are well established string- 
fragmentation models for elementary hadron hadron in- 
teractions, are included in this figure. In all these mod- 
els, the B/B ratio increases strongly with the strangeness 
content of the baryon. For strangeness \s\ = 3 the ratio 
significantly exceeds unity. In UrQMD and PYTHIA the 
hadronization of the diquark quark strings leads directly 



1 



to the overpopulation of f2. In NeXuS, however, the im- 
balance of quarks and anti quarks in the initial state leads 
to the formation of q va \ — s sea strings (the s va i — 9 sea string 
is not possible) . These strings then result in the overpop- 
ulation of fi's. 
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FIG. 1. Left: anti baryon to baryon ratio at \y — 2/ cm | < 1 in 
pp interactions at 160 GeV as given by PYTHIA, NeXuS and 
UrQMD. Right: anti baryon to baryon ratio for the same re- 
action as given by statistical models. Stars depict preliminary 
NA49 data for the B/B ratio at midrapidity. 

In Fig. |l| the string model results are compared with 
the predictions of statistical models (SM). Within the 
SM, hadron production is commonly described using 
the grand canonical (GC) partition function, where the 
charge conservation is controlled by the related chemical 
potential. In this description a net value of a given U(l) 
charge is conserved on average. However, in the limit of 
small particle multiplicities, conservation laws must be 
implemented exactly. This is done by using the canon- 
ical (C) ensemble 0,f^,f||. The conservation of quan- 
tum numbers in the canonical approach severely reduces 
the phase space available for particle production. Thus, 
exact charge conservation is of crucial importance in the 
description of particle yields in proton-induced processes 
and in e + e~ |llf| , as well as in peripheral heavy-ion col- 
lisions m. 

In Fig. [l]_the predictions of two different statistical 
models for B/B ratios in pp collisions are included. The 
main difference between these models are the implemen- 
tation of baryon number and electric charge conservation 
and the way an additional strangeness suppression is in- 
troduced. 

(I) The calculation in this statistical model |22| is a 
full canonical one with fixed baryon number, strangeness 
and electric charge identical to those of the initial state. 
Also, an extra strangeness suppression is introduced to 
reproduce the experimental multiplicities. This is done 
by considering the number of newly produced (ss) pairs 



as an additional charge to be found in the final hadrons. 
The ss pairs fluctuate according to a Poisson distribution 
and its mean number is considered as a free parameter 
to be fitted |p2| . The parameters used for the prediction 
of the f2 + /£l~ ratio (T, global volume V sum of single 
cluster volumes and (ss)) have been obtained by a fit to 
preliminary NA49 pp data Q yielding T = 183.7 ± 6.7 
MeV, VT 3 = 6.49 ± 1.33 and (ss) = 0.405 ± 0.026 with a 
X 2 /dof = 11.7/9. It must be pointed out that the D+/S1~ 
ratio is actually independent of the (ss) parameter and 
only depends on T and V (see also Fig. |J). 

(II) Here the conservation of baryon number and elec- 
tric charge is approximated by using the GC ensemble. 
Under thermal conditions at top SPS energy this approx- 
imation leads to deviations from the exact C results in 
pp collisions by at most 20%-30% . Strangeness con- 
servation is, however, implemented on the canonical level 
following the procedure proposed in |l8| . It accounts for 
strong correlations of produced strange particles due to 
constraints imposed by the locality of the conservation 
laws. In pp collisions strangeness is assumed not to be 
distributed in the whole volume of the fireball but to be 
locally strongly correlated. A correlation volume param- 
eter Vq — 47t_Rq/3 is introduced, where Ro ~ 1 frn is a 
typical scale of QCD interactions. The previous analysis 
of WA97 pA data yields: Ro ~ 1.12 fm corresponding to 
Vb — 6 fm 3 . Note that, however, hidden strange parti- 
cles are not canonically suppressed in this approach. The 
analysis of experimental data in AA collisions has shown 
that T and /ib are almost entirely determined by the 
collision energy and only depend weakly on the number 
of participants ||. The 47r results of NA49 on p/n and 
^ / Apart ratios in p-p and Pb-Pb collisions coincide within 
20%-30%. In terms of the SM this can be understood if 
T and /is in p-p and Pb-Pb collisions have similar values. 
We take T ~ 158 MeV and hb - 238 MeV as obtained 
from the SM analysis of full phase-space NA49 Pb Pb 
data ||. The volume of the fireball V ~ 17 fm 3 and the 
charge chemical potential in pp was then found to repro- 
duce the average charge and baryon number in the initial 
state. 

The predictions of the statistical models are shown in 
Fig. |l| (right). In these approaches the B/B ratio is 
seen to exhibit a significantly weaker increase with the 
strangeness content of the baryon than expected in the 
string fragmentation models. For comparison, both fig- 
ures include preliminary data on the B/B ratios obtained 
at midrapidity by the NA49 Collaboration B3 ], 



Model 


n (xitr 4 ) 


n (xl0~ 4 ) 


NeXuS 


0.48 


0.79 


PYTHIA 


0.17 


0.30 


UrQMD 


0.66 


1.05 


Canonical Model I 


0.46 


0.31 


Canonical Model II 


0.41 


0.24 



TABLE I. Predictions of different models on the 4n yield 
of Q. and Q, in pp collisions at 160 GeV. 



2 



Note that the predictions of the statistical models in 
Fig. 1 refer to full phase-space particle yields whilst mea- 
surements of B/B ratios in pp collisions have been per- 
formed at midrapidity, where they are expected to be 
the largest. Therefore, sizeable deviations of the model 
results from the data seen in Fig. 1 are to be expected. 




FIG. 2. Rapidity density of and Q in pp interactions at 
160 GeV as predicted by UrQMD, NeXuS, PYTHIA. 

The rapidity dependence of the £1 and £1 yield is pre- 
sented in Fig. |^ within different string models. The re- 
sults were calculated in pp interactions at 160 GeV within 
PYTHIA, NeXuS_and UrQMD (from top to bottom). As 
can be seen, the Q/fi ratio is largest around mid rapidity. 

The fl/Cl ratio is fairly robust - different string-model 
implementations (PYTHIA, UrQMD, NeXuS) all agree 
in their predictions within ±20%. However, as shown 
in Table I the total 4ir yields of fi'S and £Ts may dif- 



fer by a factor of 4 between the different string models. 
The statistical models give in general more consistent re- 
sults, however, deviations up to 30% are not excluded. 
In string models, the particle abundances depend on the 
parameters chosen for the fragmentation scheme, while 
in statistical models they reflect the differences between 
the ensembles chosen. Thus, the absolute yields allow 
a distinction to be made between the implementations 
once experimental data become available. We will show 
now that this is a generic feature that string models give 
£!/£! > 1, whereas statistical models yield £!/£! < 1 in pp 
interactions. 
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FIG. 3. Fragmentation of a color field into quarks and 
hadrons. While in large strings £7's and Si's are produced 
in equal abundance (a), small di' quark strings suppress fl's 
at the string ends (b), sea-s quarks enhance fi's (c). 

In order to understand the large Q/f2 values predicted 
by string models one elucidates in Fig. |^ the color flux 
tube break-up mechanism. Fig. |^ shows the fragmenta- 
tion of the color field into quark-antiquark pairs, which 
then coalesce into hadrons. While in large strings fi's 
and 0,'s are produced in equal abundance (a), low-mass 
strings in UrQMD suppress $1 production at the string 
ends (b), while in NeXuS fi's are enhanced (c). Thus, the 
microscopic method of hadronization leads to a strong 
imbalance in £!/£! ratio in low-mass strings. 

The ratio depends in a strongly non-linear fashion 
on the mass of the fragmenting string. Fig. |] shows the 
£l/£7 ratio as a function of the mass of the fragmenting 
string (i.e. different beam energies in pp). One clearly 
observes a strong enhancement of £1 production at low en- 
ergies, while for large string masses the ratio approaches 
the value of £1/0 = 1 (which should be reached in the 
limit of an infinitely long color flux tube) . 

Statistical models, on the other hand, are not able to 
yield a ratio of O/O > 1. This can be easily understood 
in the GC formalism, where the B/B ratio is very sen- 
sitive to the baryon chemical potential fis- For finite 
baryon densities, the B/B ratio will always be < 1 and 
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only in the limit of [iq = may Q/£l = 1 be approached. 
These features survive in the canonical framework, where 
the GC fugacities are replaced by ratios of partition func- 
tions |jn]Jl?^ . This is shown in Fig. || (right) where the 
ratio Q/fi in pp collisions (according to the previously 
described model I) is plotted as a function of volume 
for four different temperatures. Hence, finite size cor- 
rections in the statistical model actually lead to the op- 
posite behavior |20| in the ratio of fi/Q vs. system size 
(i.e. volume replacing string mass) to that observed in 
the fragmenting color flux tube picture. 
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FIG. 4. Left: fi/fi ratio as a function of string mass. Right: 
fi/fi ratio as a function of the volume in model I. 

In conclusion, within the fragmenting color flux tube 
models we have predicted that the f2/f2 ratio is signifi- 
cantly above unity. This is in strong contrast to statis- 
tical model results, which always imply that B / B ratios 
are below or equal to unity in proton proton reactions. 
Since this observable is accessible to NA49 measurements 
at the SPS it can provide an excellent test to distinguish 
the statistical model hadronization scenario from that of 
microscopic color-flux tube dynamics. 
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